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We investigate the use of coherent optical fields as a means of dynamically controlling the res- 
onant behaviour of a variety of composite metamaterials, wherein the metamaterial structures are 
embedded in a dispersive dielectric medium. Control and switching is implemented by coherently 
driving the resonant permittivity of the embedding medium by applied optical radiation. The effect 
of embedding Split ring resonators (SRR) in a frequency- dispersive medium with Lorentz-like dis- 
persion or with dispersion engineered by electromagnetic induced transparency (EIT) , is manifested 
in the splitting of the negative permeability band, the modified (frequency-dependent) filling frac- 
tions and dissipation factors. The modified material parameters are strongly linked to the resonant 
frequencies of the medium, while for an embedding medium exhibiting EIT, also to the strength 
and detuning of the control field. The robustness of control against the deleterious influence of 
dissipation associated with the metallic structures as well as the inhomogeneous broadening due to 
structural imperfections is demonstrated. Studies on plasmonic metamaterials that consist of metal- 
lic nanorods arranged in loops and exhibit a collective magnetic response at optical frequencies are 
presented. Control and switching in this class of plasmonic nanorod metamaterials is shown to be 
possible, for example, by embedding these arrays in a Raman active liquid like CS2 and utilizing 
the Inverse Raman Effect. 



PACS numbers: 81.05.Xj, 78.67.Pt, 41.20. Jb, 42.50. Ar 

In recent times, several new developments have consid- 
erably altered our perception of the interaction of electro- 
magnetic waves with material media in general and light- 
matter interaction in particular. On one hand, the co- 
herent control of quantum systems utilizing the quantum 
interference route [l[ has enabled researchers to demon- 
strate superluminal Q and subluminal [|| propagation 
of light, stopped light [H, Hj], render an opaque medium 
transparent at a designated frequency (Electromagneti- 
cally Induced Transparency or EIT) |(| 0] and generate 
an ultra- high index of refraction Q, while on the other, 
the development of structured composite media known 
as metamaterials @ has opened up new possibilities like 
negative refraction [l(], [H| , perfect lenses [l2[ , " invisi- 
bility cloaks" [13j ]. etc. These effects have been realised 
using on two different approaches to the study of light- 
matter interaction. While coherent control is a purely 
quantum mechanical phenomenon, metamaterials can be 
analysed within the framework of classical electromag- 
netic theory. 

One of the most attractive features of metamaterials, 
(whose resonant response is governed by their underly- 
ing geometric structure) is the fact that they can be de- 
signed to operate at any predetermined frequency. In 
addition, it is also possible to modify and control their 
response. However, the very resonant nature of the re- 
sponse of metamaterials that gives rise to their novel 
properties also hinders their applicability owing to the 
high levels of dissipation associated with their response. 
Although designing metamaterials that function at op- 
tical frequencies is difficult, it has been demonstrated 
by several groups [l4j]. Another drawback of metama- 
terials arises from the fact that since their response de- 
pends on their geometrical structure, their properties are 
fixed once they have been fabricated. However, to incor- 



porate greater flexibility to their structure and perfor- 
mance, metamaterials need to be both reconfigurable as 
well as controllable. Recently, control of metamaterial 
response has been demonstrated, both by passive as well 
as dynamical methods [TB43- 

In this paper, we demonstrate the possibility of con- 
trolling the magnetic response of metamaterials consist- 
ing of Split Ring Resonators (SRR) [l9| by actively tun- 
ing the capacitance of these structures at near-infra- 
red and optical frequencies. We have used the term 
metamaterial to signify the bare metamaterial while 
the term 'composite metamaterials' indicates metama- 
terials embedded in a dispersive background. Our ap- 
proach to this scheme is both simple and straightfor- 
ward. It makes use of quantum mechanical phenomena 
to paramctrically control a purely classical effect. This is 
achieved by embedding the metamaterial within a coher- 
ent atomic/molecular medium. An appropriate choice of 
the resonance frequencies of the metamaterial and the 
embedding medium, together with the control field and 
its detuning (in the case of EIT) results in complete con- 
trol of the effective magnetic response of the metamate- 
rial. This occurs through the interaction of the narrow 
atomic resonance of the background dielectric with the 
broader magnetic resonance of the metamaterial. Along 
with the control over the magnetic response, it is seen 
that the dissipation inherent in the metamaterial can be 
reduced significantly. This scheme may be used to pro- 
duce metamaterials which are reconfigurable as well as 
controllable by choosing a slightly different resonance fre- 
quency for the embedding medium. In addition, we also 
show how an array of nanorods that exhibits a magnetic 
response which is plasmonic in origin [20| can also be 
controlled via the self-capacitance of the array. Further, 
the same array of nanorods (which is known to behave 
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FIG. 1: Atomic level schemes for realizing (a) Lorentz per- 
mittivity dispersion(by means of a Raman transition) and (b) 
EIT-like permittivity dispersion. 

(c) The homogenizable SRR used in our calculations. A 
slab consisting of these SRRs exhibits negative magnetic per- 
meability in a frequency range above the magnetic resonance 
frequency u m = LC. The following structural parame- 
ters are used in our calculations: a = 600nm, b = 312nm, L 
= 144nm, d = 24nm, D = 24nm. 



where = Ne 2 /meo (N is the density of (atomic) os- 
cillators, eo is the permittivity of free space while e and 
m are the electronic charge and mass, respectively), u! e 
is the dielectric resonance frequency and j e is the dissi- 
pation factor. Such a dispersion can be realized by using 
for example, resonant quantum dots or choosing a Ra- 
man transition [22] with a strong pump field ensuring 
that the resonant Raman probe frequency is close to the 
magnetic resonance frequency of the metamaterial, u> m . 

In the second case, when the medium exhibits EIT, the 
permittivity for the probe field is given by [23[ 
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like an effective plasma when radiation whose electric 
field is axially oriented is incident upon it), can be made 
transmittive at frequencies well below its effective plasma 
frequency via the resonantly enhanced permittivity of its 
background medium. 



I. ANALYTICAL DESCRIPTION 

Let us consider a metamaterial consisting of an ar- 
ray of symmetric split ring resonators (SRR), as shown 
in Fig. Ufa). These arrays lie in the horizontal plane and 
are stacked up vertically. The size of the individual in- 
clusions (~ A/10) makes it possible to describe the ar- 
ray as an effective medium at the frequencies of interest. 
This medium has been shown to possess a resonant mag- 
netic response and its effective permeability is /ibaro given 
by M- 



Mbare(w) = 1 + 
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Here w\ = El h Es , fi c is the Rabi frequency of the control 

field defined as tt c = k = (N a \(d 13 )\ 2 )/(e h) is the 

strength of the transition, N a is the atomic density, 713 
and d 13 represent the decay rate and the dipole moment, 
respectively, between the atomic levels |1) and |3) while 
di2 represents the dipole moment between levels |1) and 
1 2) (see Fig. [Ijb)). It should be noted that in the absence 
of the control field, the dispersion obtained is the same 
as in the case of the Lorentz permittivity. 

In the presence of the embedding medium, the ca- 
pacitance of individual SRR units becomes strongly 
frequency-dependent. This makes the resonance fre- 
quency of the composite medium (u e fi) appear frequency 
dispersive as: ui c b — f m ■ 



The sharply resonant nature of e(w) results in the SRR 
medium showing a resonant magnetic response at mul- 
tiple frequencies. These predominantly occur above and 
below uj m , depending on whether the frequency of the 
incident wave is greater than or lesser than uj e . The new 
resonance frequencies are the solutions of the transcen- 
dental equation 



Here, io m = ^ c where L and C are the total induc- 
tance and the total capacitance of the structure [2l[ , f m 
is the filling fraction while the factor r m , which governs 
the effective dissipation present in the SRR medium, is 
related to the intrinsic dissipation rate of the metal de- 
fined according to the Drude model for the dielectric per- 
mittivity. This medium exhibits a negative permeability 
leading to a band gap in a frequency range above u m if 
r m is small enough. We study the effect of embedding 
an array of such SRRs in a medium having a frequency- 
dispersive dielectric permittivity. This may either be a 
resonant Lorentz-type dispersion, or under the applica- 
tion of a pair of coherent control and probe fields, EIT- 
like. In the first case, the relative permittivity of the 
medium is given by the expression 
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while Im[wJ,/e(w)] quantifies the dissipation in the effec- 
tive medium. Essentially, the presence of the dispersive 
dielectric permittivity of the embedding medium turns 
the SRR into a resonantly driven, actively tuned LC- 
circuit. The capacitance of the of the LC-circuit becomes 
frequency dependent, allowing the satisfaction of the res- 
onance conditions Eq. Q at multiple frequencies with 
the corresponding negative-/! band gap of the 'bare' SRR 
splitting into two or more band gaps. 

The filling fraction f e g and the dissipation parame- 
ter of the composite metamaterial can be viewed as ef- 
fective frequency dependent quantities whose behaviours 
are greatly influenced by the resonant background. The 
effective filling fraction f e g and the effective dissipation 
parameter T e g have the following generic forms in the 
presence of a frequency-dispersive background: 
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FIG. 2: (a)The frequency dependent (j, e g and (b) r o fr and 
filling fraction f e a (obtained analytically) for the SRR meta- 
material in the presence of a embedding dielectric medium 
displaying EIT. (c) and (d): The same quantities for the SRR 
metamaterial embedded in a medium which has a resonant 
Lorentz permittivity. 



current induced in the SRR unit by the magnetic field of 
the incident radiation. The second arises due to electric 
polarization of the embedding dielectric background by 
the electric field of the incident light. At resonance, the 
two are out of phase with each other leading to the freez- 
ing of currents in the SRR. The proximity of u e to co m 
determines the extent of modulation of the )i e e response. 
In other words, the interaction between the electric and 
magnetic resonances is the strongest when u e lies in the 
vicinity of u) m , (w m ±r m ). Otherwise, the two resonances 
are virtually uncoupled. 

The magnetic response of the composite metamate- 
rial embedded in an EIT medium exhibits multiple mag- 
netic resonances (Fig. [UJb)). This is again accompanied 
by drastically reduced absorption within two frequency 
bands. The strength and the detuning of the control field 
can both be used to shift the positions, widths, slopes 
and dissipation associated with the various resonances in 
/i c ff. In fact, this scheme is more amenable to control as 
compared to the one described earlier because the detun- 
ings of the probe and the control beams (both positive 
and negative) can be used to shift the EIT line centre 
and hence, the frequencies at which the absorption in 
the metamaterial is reduced. Again, the effective filling 
fraction and the modified dissipation parameter are given 
by 



/off M = fm [{Real [e(u)}} 2 + {Im [e(oj)}} 2 ] , (5) 

r cff H = ^J&A + r m {Real [eH]} 2 

+r m {Im[eH]} 2 . (6) 

Let us now study the effect of the Lorentz and the 
EIT dispersion on the properties of the metamaterial in 
some detail. If the background dielectric medium has a 
Lorentz-type dispersion, the following expressions for f e g 
and r e ff are obtained: 
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Along with the two new magnetic resonance frequencies 
and a split in the negative /Lt-band, a region of signifi- 
cantly reduced absorption is observed between the two 
resonances. In fact, the Im[// e g] is observed to pass 
through a minimum at the resonant frequency of the bare 
SRR if co m = io e . There is also a significant change in the 
nature of f e s and T e g around the frequency u e . The 
large values of f e g at ui e is due to the resonant enhance- 
ment of the cross-section of the SRR structural unit: the 
SRR appear much larger to radiation at resonance than 
they actually are. The reduction in absorption at ui e is 
due to the vanishing of the currents in the SRR. This 
results from an interplay of two currents giving rise to 
different charge densities. The first of these is due to the 
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,(10) 



where A = lj\ — ui. 

From the preceding expressions, we find that the new 
filling fraction f e g increases close to the resonant fre- 
quency while r e ff exhibits sub-natural (r c ff <C T m ) val- 
ues around the atomic resonances (Fig. [He)). At the 
EIT line centre frequency (wi), f c g and T c g are found 
to have the same values as the bare SRR medium. As 
before, the reduction in absorption is a manifestation of 
the vanishing currents in the SRR loops. It is interesting 
to note that if losses are reduced significantly, the regions 
where Re [/x e ff] < gradually disappear, (see Fig. [2]), in 
agreement with the findings of Stockman [24j . 

If the metamaterial has been designed to exhibit a res- 
onant response at the frequency of the EIT line centre in 
the absence of medium, the resonant magnetic response 
of the composite metamaterial is found to be such that 
Re [/icff] — > 1 (with an associated imaginary part) ex- 
actly at the EIT line centre. However, if uj m ^ the 
response of the composite at the EIT line centre is simply 
the value of ^(bare at the frequency uii . It should be em- 
phasized that since we are dealing with a resonant system 
with a highly dispersive lineshape, changes of the control 
parameter (ej, or ceit) can lead to large-scale changes in 
response (/i c ff) at a given frequency, due to the shift in 
the resonance frequency. 
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FIG. 3: (a): The computed band structures for the bare SRR 
(red circles) and the composite metamaterial embedded in a 
Raman medium (black squares) under the condition uj e = tj m . 
The blue lines represent the dispersion predicted by the an- 
alytic formula. The band gaps due to (j, e g < are indicated 
by the cross-hatched region on the right (for the bare SRR 
whose cum = 74.9THz) and the the two hatched regions on 
the the left (for the band gaps resulting from the inclusion of 
the resonant Raman medium). 

(b) The reflectivity and the transmittivity for a calculated for 
a slab consisting of a single layer of SRR. Note the presence 
of two stop bands (low \T\ 2 ) corresponding to the band gap 
frequencies in (a) and the new transmittive band that devel- 
ops within the original band gap. 

(c) The effective impedance Z(//)of the composite metamate- 
rial showing the presence of two resonances. 

(d) : The /i e ff, for a slab of the SRR metamaterial with an 
embedded resonant medium, retrieved using the computed re- 
flection and the transmission coefficients. Note the existence 
of two frequency ranges where fj, e g < 0. The grey dotted line 
represents the zero level. 



II. NUMERICAL SIMULATION OF THE SRR 
METAMATERIAL AT MID-IR FREQUENCIES 

Numerical simulations using the PHOTON codes 
based on the Transfer Matrix Method [25|, [26| were car- 
ried out to verify our theoretical predictions. We present 
the canonical example of subwavelength-sized SRR made 
of silver at mid-infrared (mid-IR) frequencies. The meta- 
material shows magnetic activity for applied radiation 
polarized with the magnetic field along the SRR axis. 
The structures were modelled as being embedded in a 
resonant background. We have used experimentally ob- 
tained parameters for silver [27] in our calculations. (We 
use silver for all metallic inclusions in this paper.) 

The cross-section of the (cylindrical) SRR that we con- 
sidered in this case is shown in Fig. [T] Seventy-five grid 
points along each direction were used for computing the 
band structure as well as the reflection and the transmis- 
sion coefficients. Invariance along the cylindrical axis is 
assumed. The bare SRR in this case has a magnetic reso- 
nance frequency of 74.9 THz. This metamaterial contain- 



ing a Lorentz-type dielectric within its capacitive gaps is 
modelled using the following parameters: f e — 24.1 THz, 
iu e = 74.9 THz and 7 e = 2.4THz. 

As predicted in the preceding section, the negative fi 
band splits in the presence of the embedding dielectric 
medium. The frequencies at which the new band gaps 
are formed are consistent with the predicted values. The 
mismatch between the results of the TMM calculation 
and the analytic formula can be attributed to the pres- 
ence of parasitic capacitances in the SRR. These capac- 
itances, which were neglected in the theoretical model 
used to describe the functioning of the SRR [2l| , are au- 
tomatically taken into account in the TMM simulation. 

The results of the band structure calculation were com- 
pared with the transmittivity and the reflectivity calcu- 
lated for a slab consisting of a single layer of the SRR. 
A new propagating band (indicated by a peak in the 
transmittivity) is found to occur within the stop band 
region of the bare SRR. The frequency intervals corre- 
sponding to the new band gaps (regions of high reflectiv- 
ity and low transmittivity) were found to be consistent 
with the band structure calculations. To illustrate the 
splitting of the negative permeability band, the effective 
material parameters for a slab consisting of four layers 
of the SRR-dielectric composite were determined using 
a retrieval procedure that utilizes the complex reflection 
and transmission coefficients (28|. The splitting of the 
negative /i band is found to be consistent with the pre- 
dicted values (Fig. [3] (d)) and also with the simulated 
band structure and the reflectivity and the transmittiv- 
ity. In addition, Fig. [3f c) shows the effective impedance 
Z(i/) of the metamaterial, with two resonances at the 
same frequencies as the negative permeability bands, in 
contrast with the single peak in 1(y) observed for the 
bare SRR[U . 

Fig. @] illustrates the response of the metamaterial 
when the dielectric and the magnetic resonances do not 
coincide, i.e., uj e ^ co m . As before, the negative per- 
meability band splits, but the new resonances are either 
red-shifted or blue-shifted (with respect to those obtained 
when u e = u) m ), depending on whether ui e < uj m or 
uj e > uj m . A study of the reflection and transmission 
coefficients for a slab consisting of four layers of the SRR 
further confirms this behaviour. As seen in the case when 
uj m = w e , two regions of enhanced reflectivity (corre- 
sponding to the band gaps) separated by a transmittive 
region (corresponding to the additional propagating band 
that develops) are observed. For multiple layers, addi- 
tional peaks due to Fabry-Perot resonances (arising as a 
result of multiple scattering) are observed. In each case, 
the calculated effective impedance (Z(i/)) of the structure 
(not shown here) reveals the presence of two peaks, sig- 
nifying the two magnetic resonances. However, the effect 
of the resonant Lorentz medium weakens considerably 
when the two resonances are far apart. 

In a similar fashion, appropriately designed SRR struc- 
tures, embedded in resonant atomic media, can be 
shown to behave in the same way at NIR/optical fre- 
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FIG. 4: The band structures and (|T| 2 , |_R| 2 ) for the compos- 
ite metamaterial when u m 7^ ui e . 

(a) : Band structure of the composite metamaterial when w e 
= 79.69THz. The resulting gaps are blue-shifted with respect 
to those obtained when uj e = ui m . 

(b) : Band structure of the composite metamaterial when lj £ 
= 70.03THz. In this case, the resulting gaps are red-shifted 
with respect to those obtained when cu e = cj m . 

(c) & (d): Reflectivity and transmittivity calculated for a slab 
consisting minima in \T\ 2 correspond to the new band gaps. 



quencies [29(. Multiple layers of such SRR can switch 
sharply between highly reflecting to transmitting states 
within a narrow bandwidth. These fine features in the 
|-R| 2 — \T\ 2 profile arise due to the Fabry-Perot resonances 
of the slab. Coupled with the control possible via EIT, 
these materials are potential candidates for narrow-band 
switching applications. 



III. ROBUSTNESS OF THE CONTROL 

Although it is often seen that light interacting with 
mesoscopic metallic structures experience much higher 
dissipation than the bulk metal itself, we have found that 
our results remain unaffected even if the intrinsic dissi- 
pation of the metal (silver, in this case) is increased three 
times from its bulk value. 

There is another aspect of these structures that needs 
to be taken into account. The resonance frequency of 
these metamaterials depends crucially on the accuracy 
with which they can be fabricated. Thus, in prac- 
tice, each metamaterial unit will be different from the 
rest, leading to a distribution of resonance frequencies 
rather than a single resonance frequency. The consequent 
broadening of the magnetic resonance could wipe out the 
narrow band effects predicted. Fortunately, we find that 
the effects are rather robust against the such inhomoge- 
neous broadening. 

We consider a distribution of magnetic resonance fre- 
quencies arising due to structural imperfections to be nor- 
mally distributed (Gaussian) about the mean frequency 
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FIG. 5: (a)The effective permeability of the metamaterial for 
different widths of the distribution of the magnetic resonance 
frequency. The resonant fi e g is washed out when the imper- 
fections lead to a ~ u m . 

(b):The broadened Im(^ e ff) of the metamaterial in the pres- 
ence of a background medium which exhibits EIT, for a = 
cj m /I00. The inset shows the details of the central part of 
the figure where the narrow features are seen to survive in- 
spite of the broadening of the overall response. 



oj m . The width of this distribution (er) is indicative of 
the extent of the imperfect structuring. Averaging over 
this distribution, we obtain: 



(/off') 



^eff 

r c ff 

/cff 



(11) 

(12) 
(13) 



Note that the effective filling fraction (/ e ff) remains un- 
altered. Thus, the effective permeability of the medium 
is found to broaden, as shown in Fig. EJa). 

For values of a ~ T m , we find that the narrow band ef- 
fects arising due to the presence of e(u>) survive (See Fig. 
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[5]Jb)). While the strength of the resonance decreases on 
the whole due to the larger V e g, the regions of reduced 
absorption between the resonances remain nearly unal- 
tered. The 'increased' filling fraction that remains un- 
affected by the inhomogeneous broadening compensates 
for the effect of the imperfections. In other words, the 
few resonant SRR units exhibit a large enough absorption 
cross-section such that the reduced number of participat- 
ing SRRs due to the inhomogeneous broadening does not 
adversely affect the response on an average. Although 
accuracy in fabrication still remains an important issue, 
we note that the controllable effects discussed here are 
robust and generic to many metamaterials that depend 
upon resonances for their electromagnetic properties. 



IV. PLASMONIC METAMATERIAL AT 
VISIBLE FREQUENCIES 

In this section, we present the response and control 
of a metamaterial consisting of plasmonic nanorod loop 
inclusions. Arrays of small resonant plasmonic spheres 
arranged on subwavelength-sized loops have been shown 
to exhibit a magnetic response in Ref. [20| . We consider 
here plasmonic inclusions of metallic (silver) nanorods 
arranged on the circumference of a circle (Fig. (Hla)). 
The magnetic field of the incident radiation is assumed 
to be oriented along the axis of the rods. On assum- 
ing invariance along the rod axis, the problem becomes 
essentially two-dimensional. The loop exhibits a mag- 
netic resonance that stems from an interaction driven 
by the individual plasmonic resonance of each nanorod. 
Each resonating nanoparticle induces a circulating dis- 
placement current around this loop. The electric dipole 
moments (perpendicular to the axis) tend to align along 
the circumference of the ring when driven at the plasmon 
resonance frequency. This is confirmed by our COMSOL 
finite elements numerical simulations (Fig. [Sib)) where 
an enormous confined magnetic field appears in the inte- 
rior of the loop at resonance, when driven by an external 
source. 

The size of the loop does not strongly influence the 
frequency of the individual magnetic resonance occurs. 
The plasmonic resonance frequency of each nanoparticle 
primarily determines the resonant behaviour of the loop. 
However, the number of plasmonic particles on the loop 
determines the strength of the resonance [2(|. At res- 
onance, the size of each unit cell (a ~ A/5) and each 
loop inclusion is much smaller compared to the wave- 
length of the incident light (i? ~ A/25). Thus, Bragg 
scattering does not play a dominant role in the observed 
phenomenon and the system is easily homogenized. 

The response of an array of such plasmonic loops was 
calculated using the PHOTON codes. The calculations 
are essentially two dimensional, the cylindrical axis being 
the axis of invariance. To avoid staircasing effects with a 
square grid, we considered rods with square cross-sections 
rather than circular ones. The unit cell is a square with 




FIG. 6: (a): A cross-sectional view of a typical unit cell 
of the plasmonic loop metamaterial. The picture shows the 
nanorods placed on the circumference of a circle. The grey 
area represents the embedding dielectric medium, 
(b): The magnetic response of an isolated loop inclusion il- 
lustrating the concentration of the magnetic field inside the 
loop at resonance, when excited by means of a line source. 
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FIG. 7: (a) The band structure and (b), the transmittivity for 
four layers of unit cells of the plasmonic metamaterial. The 
additional peaks are the Fabry-Perot resonances occuring due 
to multiple scattering between the interface of the medium 
and vacuum. 

side a = I50nm (See Fig. [U(a)). Within each unit cell of 
the metamaterial, four silver nanorods (of square cross- 
section and side b = 25nm ), embedded in a background 
dielectric medium with eb = 2.65, are symmetrically ar- 
ranged on the circumference of a circle of radius R = 26 
nm. The metamaterial exhibits a distinctly identifiable 
magnetic resonance at ~ 415 THz, with a negative mag- 
netic permeability band gap up to ~ 421 THz. This be- 
haviour is also confirmed by the transmittivity for a slab 
composed of four layers of unit cells (loops) (Fig. [7f b)). 
We note that the resonance frequency of the metama- 
terial differs from that of the individual loop inclusions 
due to coupling between the nanoparticles. For smaller 
unit cell sizes and larger loop radii, the coupling between 
neighbouring plasmonic loop resonances leads to the hy- 
bridization of bands. Many weak resonances are found 
to occur close together and identification of the magnetic 
resonance leading to the negative permeability band gap, 
is complex for larger loop radii. 

It is well-known that the plasmonic resonance of a 
metallic nanoparticle is highly sensitive to the dielec- 
tric permittivity of the surrounding media (e;,). The 
plasmonic resonance of a cylindrical nanoparticle occurs 
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FIG. 8: The proposed level scheme for Inverse Raman Effect 
in CS 2 . 
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FIG. 9: (a) & (b) The new band formed and the modified re- 
flectivity of the plasmonic metamaterial in the neighbourhood 
of a Raman transition in CS2 at 416.54 THz. The vertical line 
in (b) represents the reflectivity of the bare 'metamaterial'. 



when the condition 



Re[e b 



J=0, 



(14) 



is satisfied [3C|. Using a Drude form for the permittivity 
of the metal (e m ), the plasmon resonance frequency of a 
metallic cylinder may be written as: 



1 



7 



1/2 



(15) 



where u r is the plasmon resonance frequency for the 
nanoparticle, w p and 7 defines the bulk plasma frequency 
and the bulk dissipation of the metal, while Coo is its 
'static' permittivity. Hence, even a small change in ef, 
can shift the plasmon resonance frequency of the nanorod 
and by extension, the magnetic loop resonance frequency 
arising out of the circulating displacement current around 
the loop. We note a subtle difference between the SRR 
and the plasmonic metamaterial which may affect the 
efficiency with which their magnetic response is manipu- 
lated. In the case of the SRR, the displacement current 
is confined to the region within the capacitive gaps of the 
SRR. In this case, however, the displacement current is 
between adjacent nanoparticles, and hence, spatially de- 
localized. A more closely packed loop would, in principle, 
enable a higher degree of localization of currents and a 
stronger resonant response. 

Control of e& by an applied electromagnetic radiation 
is the key to the control of the composite metamaterial 
itself. In order to demonstrate parametric control over 
the magnetic resonance, an atomic/molecular system ex- 
hibiting a resonance in the neighbourhood of the band 
gap frequency is essential. Raman transitions that pro- 
vide such resonances can be effectively used in a very 
flexible manner because the Raman resonance frequency 
can be tuned by the choice of the pump laser frequency. 
A strong pump field with frequency vq will create con- 
ditions for the resonant absorption of a probe field at 
a frequency (v + vr) (Inverse Raman Effect [3l], HH), 
where vr is the molecular level spacing of the medium. 
We present the calculated results for such a stimulated 



Raman absorption process (Fig. [8]) at 416.54 THz in CS2 
with ground state vibrational level spacing of 656 cm -1 
(19.68 THz). Experimentally obtained values for the pa- 
rameters of CS2 were used [33j . corresponding to the 
non-resonant response. We use peak pump pow- 
ers of the order of 100 GW/cm 2 and a typical linewidth 
of 0.5cm _1 (15 GHz). Our estimates for the nonlinearity 
are overtly conservative as the resonant Raman process 
would be expected to enhance the x response. Such 
inverse Raman absorption of about -15dB has been re- 
ported in amorphous silicon by with much lower pump 
powers of about 1 GW/cm 2 [351 ] . assuming a pump laser 
field at 396.86 THz which lies within the propagating 
band of the metamaterial (Fig. [TJa)). 

An additional propagating band develops within the 
band gap region seen in Fig. in the presence of the 
pump field and this propagating band region has been 
expanded in Fig. |H] This is manifested by a sharp drop 
of (~ 20%) in the reflectivity of the medium. The new 
band formed is quite dissipative and this prevents a 
large change in the transmittivity in comparison with 
the change in reflectivity. The transmittivity of the com- 
posite metamaterial is reasonable when compared to that 
for the 'bare' metamaterial, but is completely overshad- 
owed by the large change in the reflectivity. It is pos- 
sible to numerically demonstrate dramatic enhancement 
of the transmittivity by several orders of magnitude by 
using higher pump intensities. Thus, the effect of the 
introduction of the resonant Raman medium results in a 
sharp reduction in the reflectivity of the composite meta- 
material over a narrow frequency band (~ 5 GHz). It is 
the linewidth of the Raman pump laser that dictates the 
linewidth of the induced absorption, which in turn deter- 
mines the change in the reflectivity. The number of prop- 
agating bands that are formed are strongly dependent on 
the resonant nature of the background medium. In the 
case of EIT-based control, additional bands which switch 
rapidly between states of low and high transmittivity, are 
formed in the presence of an control field fi c [2^] . 

The behaviour of the metamaterial for radiation of 
the other state of polarization (with the electric field 
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of the incident radiation directed along the axis of the 
metallic rods) may also be controlled in a similar man- 
ner [34j]. Fig.|6](a) continues to represent a typical section 
of the metamaterial, considered here, however, with the 
E field along the nanorod axes. As shown in Ref. p34j . 
the transmittivity of the effective plasma is zero while 
the reflectivity is high below the cut-off frequency. At 
higher frequencies, the transmittivity of the medium in- 
creases considerably. These correspond to the regions 
where e c g < and e c g > 0, respectively. The response of 
this metamaterial can be actively manipulated using em- 
bedded atomic/molecular media. The introduction of an 
embedding medium with a resonant response below the 
effective plasma frequency, into the metamaterial, will 
drive the effective permittivity to positive values over 
narrow frequency bands. This results in the formation 
of transmittive bands below u> p . Resonant enhancement 
of the background permittivity can also be achieved by 
stimulated Raman absorption due to an applied pump 
radiation in a similar manner. A new transmission band 
is found to develop at this frequency, which lies below the 
effective plasma frequency for the medium [34]]. We fur- 
ther note that a lower linewidth of the pump laser results 
in a significant increase in transmittivity. The plasma- 
like behaviour of the nanorod array depends on a bulk 
averaging process to determine its overall response and 
does not necessitate a 'loop'-like arrangement. Thus, the 
response of the nanorod array to light of the two orthog- 



onal states of linear polarization and their subsequent 
modified behaviour in the presence of a dispersive per- 
mittivity has a different origin. 



V. CONCLUSION 

We have thus demonstrated a new scheme for the con- 
trol of metamaterials where the resonant metamaterial 
response is itself completely transformed by the introduc- 
tion of controllable resonant media into the design. Cou- 
pled with the fact that metamaterials can be designed to 
show a resonant response at a predetermined frequency, 
this technique provides an entirely new way of reducing 
the dissipation in metamaterials, as our results show. We 
have discussed the control of SRR-based metamaterials 
by the parametric control of the capacitance, the control 
of nanorod plasmonic loop metamaterials by the control 
of the surface plasmon resonance via the background di- 
electric environment and the control of the plasma-like 
behaviour of the nanorod metamaterial by changing the 
bulk average dielectric permittivity. The variety of such 
possibilities of control demonstrate the generality of this 
approach. The control scheme presented here enhances 
the tunability of metamaterials and can be effectively 
used for a variety of narrow-band applications, thus mov- 
ing towards dynamically controllable metamaterials. 
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